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Understanding high-temperature superconductivity requires a prior knowledge of the 
nature of the enigmatic pseudogap metallic state |jl||, out of which the superconduct¬ 
ing state condenses. In addition to the electronic orders involving charge degrees of 


a magnetic intra-unit-cell 
to set in just at the pseu- 


freedoms recently reported inside the pseudogap stat e Il2| 

(lUC) order was discovered in various cuprates jll, 0, 0, 18 
dogap temperature, T*. In nearly optimally doped YBa2Cu306.85? polarized neutron 
scattering measurements, carried out on two different spectrometers, reveal new fea¬ 
tures. The order is made of finite size planar domains, hardly correlated along the 
c-axis. At high temperature, only the out-of-plane magnetic components correlate, 
revealing a strong Ising anistropy, as originally predicted in the loop current model 


Below T*, a correlated in-plane response develops, giving rise the apparent tilt 


of the magnetic moment at low temperature fl], 0]. The discovery of these two regimes 
put stringent constraints on the intrinsict nature of lUC order, tightly bound to the 
pseudogap physics. 

Polarized neutron scattering e^eriments reported the existence of an lUG magnetic order in four 
different cuprates families [l|, [j, y-l9j . The lUG magnetic order that breaks time-reversal symmetry 
(TRS) but preserves the lattice translational symmetry (LTS). The lUG order develops below a tem- 
pCTature T^ag that matches the pseudogap temperature T* as dehned by the resistivity measurement 
|l4j. IlSj l and where recent resonant ultrasound spectroscopy measurements in YBa2Gu306+x show 
that the pseudogap state is a true symmetry breaking phase [l^. While bulk neutron scattering 


measurements provides unambigous evidence for an lUG magnetic order, such an order was not 
detected by local magnetic probe measurements jsi [13, 18|, a conundrum that could be understood 
in terms of hnite size and slowly fluctuating magnetic domains 19|], whose existence has not been 
established so far. 

The observed lUG magnetism can simply be described as a q=0 antiferromagnetic state as it 
keeps the LTS, but requires an internal staggered magnetic pattern inside each unit cell. Such a 
state is qualitatively consistent with the loop current (LG) model for the pseudo-gap [13, 19|. While 
out-of-plane magnetic moments should be produced by the planar conhnement of LG, the observed 
magnetic moments also display an unexpected in-plane component, giving rise to a tilt of about ~ 45° 
with respect to c-axis [l|, l2|, a, Isj at low temperature. Several attempts to overcome this discrepancy 
were proposed: an additional spin response through spin-orbit coupling Bol. a geometrical tilt of 
moment through delocalisation of the LG over GuO octaedra or pyramids (2l|, l22|, a superposition of 
LG states through quantum effects |23|. Besides the LG model, electronic pha ses breaking TRS, but 
preserving LTS, were proposed to preempt either a pair density wave state [2^ or a composite charge 
density wave state |25|. Alternative scenarii consider an lUG order based on a nematic-like state 
with either spin [l| or orbital j26| moments on oxygens, or the entangled spin and orbital degrees 
of freedom within a model of magneto-electric multipoles [1^. All of these theories give different 
interpretations of the origin of the tilt of the magnetic moment. 
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A this stage, more experimental information are reqnired on the size and lifetime of the magnetic 
domains and exact orientation of the magnetic moment to make progress on the nnderstanding of 
the intrinc natnre of the lUC magnetic order and its interplay with the psendo-gap physics. This 
can be achieved aronnd optimal doping. Indeed, a fast decay of the magnetic intensity was re- 
orted from the underdoped to optimally doped samples observed in Bi 2 Sr 2 CaCu 2 08+5 compounds 
Meanwhile, the magnetic intensity is proportional to Tmag in the underdoped regime, whereas 
this scaling suddenly breaks down upon approaching optimal doping. These effects could actually 
be triggered by a redistribution of the magnetic scattering in momentum space, namely a short¬ 
ening of the magnetic correlation length upon increasing hole doping. In nearly optimally doped 
YBa2Cu306.85 (Tc = 89K, p = 0.15), these characteristic properties are present (see supplementary 
materials): the lUC magnetic order still develops at a rather high temperature T^ag ~ 200K, but 
the measured intensity at the Bragg position is strongly reduced, as compared to more underdoped 
samples. 

To gain a new insight of the lUC magnetic order, we carried out on YBa2Cu306,85 a polarized 
neutron experiment using the multi-detector diffractometer D7, usually used to measure diffuse 
magnetic scattering. A survey of a few trajectories around different (1,0,L) Q-positions reveals the 
persistence of magnetic scattering away from the Bragg reflections. In particular, for the trajectory 
across Q = (1,0,0.25), the contrast in intensities noticeably improves compared to the Bragg position 
Q = (1,0,0). In Fig. [Ha-b, measurements around Q = (1,0,0.25) show that the magnetic signal 
increases upon cooling. At high temperature (300K) the magnetic signal is assumed to be flat 
whereas at low temperature (lOOK) there is a net enhancement of the magnetic intensity. From 
the H-dependences of Fig. [TJb, one can extract the q-width (full width at half maximum) by a 
Gaussian htting, Aq = 0.025 ± 0.004 r.l.u. We observe no noticeable evolution of that q-width 
with temperature. Actually, Aq is slightly broader than the resolution Aq^es = 0.0195 ± 0.001 r.l.u. 
suggesting short-range correlations. One can extract a finite planar magnetic correlation length, 
^ab = a/(7r * Aqi) where Aqi = i/Aq^ — Aq^es^ is obtained by deconvolution from the resolution: 
^ab ~ 20a ~ 75 A (a = 3.85A is the in-plane lattice parameter). Interestingly, finite magnetic 
correlation lengths generally imply the existence of a finite characteristic time scale. This may 
explain why the observation of the lUC magnetic order has not wt been corroborated by magnetic 
local probe measurements such as nuclear magnetic resonance [l9|, due to their much longer time 
scales. In such a line of thought, the correlation length ought to stay hnite even at lower doping to 
explain the absence of a magnetic signal in local probes. Finally, it is also worth pointing out that 
X-ray measurements recently revealed charge density order (CDW) j^, 30| which appears at lower 
temperature Tcdw < T^ag with correlation lengths .^cdw 8a < .^ab 20a. That suggests the 
pre-eminence of the lUC order over the CDW instability. 


_ N 

Next, we average the intensity U of N=10 detectors {I = Y, h/N) around each Q = (1,0,L) 

i=l 

position (corresponding to different rocking sample angles) to further improve the contrast of the 
magnetic intensity. The total magnetic intensity summed over 10 detectors is reported for two L 
values not at the Bragg position as a function of temperature (Fig. [2Ja). First, it reveals that the 
magnetic correlations start actually at a high temperature clearly above T^ag ~ 200K indicative of 
pre-transitional magnetic scattering. This is already depicted in Fig. [Ha from H-scans. Second, the 
comparison of the total magnetic intensity at both L values (L=0.25 and L=0.5) in Fig. [2la shows 
that the magnetic intensities grow hrst similarly at both L and then more rapidly for L=0.25. That 
behavior suggests a change of the correlations along c* upon cooling with a gradual build-up as 
expected on approaching a phase transition. Using absolute units conversion, one can compare data 
from both 4F1 and D7 (Fig. [2lb), from which one can estimate Aq = 0.65 ±0.05 r.l.u. corresponding 
to a very short correlation length along c*, ~ 0.5c, at lOOK. It is worth to emphasize that such 
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broad L-dependence can explain the sharp decreasing of the peaked magnetic intensity versns doping, 
the L-integrated intensity being in constrast consistent with the ratio of both T^ag- 
One can also use the temperature dependence of the averaged intensity in 10 detectors to estimate 
the directional components of the magnetic moment. In Fig. |2jc-d, we report the temperature 
evolution of the the out-of-plane component Me and the in-plane component Mab- Interestingly, 
both cross-sections exhibit a noticeable distinct temperature dependences. At high temperature, the 
in-plane component vanishes, meaning that the correlated magnetic fluctuations <M^> essentially 
correspond to the out-of-plane component. That strong Ising character (along c*) remains down to 
Tmag- Below Tjnag, the in-plane component is increasing as well. As / oc lo/^ab (where Jq is the 
magnetic peak intensity), it indicates either that both magnetic components have a distinct reponse 
in temperature of their amplitude or that their correlation lengths exhibit a significantly different 
thermal evolution. Interestingly, the Ising character of the lUC order is masked at low temperature 
by the onset of the additional in-plane component. 

At lOOK, one can estimate the ratio between Mab and Me. That defines an apparent tilt angle, 6 of 
the magnetic moment with respect to the c*-axis as tan(0) = Mab/Mc jl]: from Fig. Ob-c, tan(6') ~ 
T8T JThat gives 6 ~ 40° ± 9° for L ~ 0 similar to what has been deduced previously for lower doping 
That result dismisses all models where 6 goes to zero at L=0 such as LC models that include 


apical oxygens as an explanation for the tilt [21|, [22|. The observed temperature dependence of 9 is 


also not consistent with theories where the tilt direction is due to rigid geometrical factors. 

Going further, we applied the standard D7 polarization analysis relation for paramagnetic systems, 
typically valid for disordered magnetism [lit] to get Imag around H=0.88 and L=0 as a function of 
temperature (Fig. [31a). Interestingly, one observes a maximum in this quantity at around 200K 
corresponding to the transition temperature T^ag (see Fig. S2.c-d measured on 4F1). At wave- 
vector (H,0,0) far away from H=1 (say H < 0.8), this cusp-like shape disappears. Such a shape is 
usually characteristic of critical scattering around the ordering temperature. However, here none of 
the magnetic cross-sections at H=0.88 displays any cusp as shown by the temperature dependence 
of the SF cross sections for two polarizations (Figs. [3lb). This absence is actually not surprising as 
we observe finite correlation lengths. 

This implies a more subtle interpretation of the maximum in terms of the magnetic components, 
<M^> and <M^^>, from which, one can rewrite Imag ~ 2(< M^^ > — < M^ >). At 300K, both 
components are fluctuating. When <M^> starts to correlate at H=1 (Figs. Oc), the intensity at 
H=0.88 in the Z polarization decreases accordingly whereas <M^^> remains fluctuating (Figs. Ob), 
yielding an increase of Imag- Below 200K, <M^^> correlates as well and Imag decreases accordingly. 
The cusp observation then further points toward a distinct behaviour for both magnetic components. 


Methods: 

We here report polarized neutron measurements on a YBa2Cu306.85 (Tc = 89K) single crystal, 
previously used to study the spin dynamics jl^, extending the study of the lUC order at a higher 
hole doping p = 0.15 near optimal doping. The polarized neutron experiments have been performed 
on two spectrometers: the triple axis spectrometer 4F1 (Orphee, Saclay) and the multi-detector 
diffractometer D7 (ILL, Grenoble) (see supplementary informations). For both, a polarizing super¬ 
mirror (bender) and a Mezei flipper are inserted on the incoming neutron beam in order to select 
neutrons with a given spin. In addition, a filter (pyrolytic graphite for 4F1 or beryllium for D7) is put 
before the bender to remove high harmonics. After the sample, the final polarization, P, is analyzed 
by an Heusler analyzer on 4F1 whereas D7 is equipped with an array of polarizing benders in front 
of the detectors. For each wave vector Q, the scattered intensity is measured in both spin-flip (SF) 
and non-spin-flip (NSF) channels. For the measurements on 4F1 the incident and hnal neutron wave 
vector are set to kj = kp = 2.57A“^. On the diffractometer D7, the incident wave-vector is taken as 
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ki = I. 29 A As we are using cold neutrons, we had about 4 times better Q-resolution compared 
to 4F1. Following previous studies 0,3119 , the search for magnetic order in the pseudo-gap phase 
is performed on Bragg reflections Q = (1,0,L) with integer L values. The sample was then aligned 
in the scattering plane (1,0,0)/(0,0,1). 

D7 is equipped with a hxed XYZ polarization mode[TH with Z being vertical and X and Y are 
pointing along arbitrary in-plane directions. For all of the three polarizations X, Y and Z, the SF and 
NSF channels were systematically measured over a few rocking angles of the sample in order to r nap 
the magnetic scattering. We then performed data reduction adapting the standard procedure [ul 
to force the magnetic scattering to be zero at 300K for any polarization direction. This is made in 
order to find evidence for small magnetic intensities (see supplementary information for details). The 
procedure includes background subtraction, flipping ratio, vanadium corrections. Further, following 
the procedure discussed above for 4F1, the conversion in absolute units has been performed using 
the Bragg peak Q=(1,0,0) intensity. 
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FIG. 1: (color online) In-plane magnetic correlations, (a) Three dimensional plot of the maximum of the magnetic intensity 
in the SF channel around Q = (1,0,0.25) as a function of temperature and wave-vector, (b) Cuts of the 3D-map at three 
temperatures: lOOK (red squares), 200K (brown squares) and 300K (black empty squares). Data obtained in channels X and Z 
are here averaged, c) schematic view of the reciprocal space probed on D7 with in red the measured trajectory going through 
(1,0,0.25). r.I.u. stands for reduced lattice units. 
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FIG. 2: (color online) (a) Temperature dependences of the total magnetic intensity, Ix -I-Ty +Iz (background subtracted) for 
two L. (b) Magnetic intensity in absolute units from 4F1 (empty symbols) and D7 (full symbols) versus L at lOOK. Temperature 
dependence of the magnetic components for L=0.25: (c) the out-of-plane component and (d) the in-plane component 
(see supplementary information). Solid lines are guides to the eye. 
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FIG. 3: (color online) (a) Polarization analysis around Q ~ (0.88,0,0) of Imag = (2Ix^ + 2l|^ —4l|^). The majcimum of intensity, 
which is only for Q ~ (0.88,0,0), appears around Tmag 200K. (b) Temperatnre dependences of the SF cross section in the 
(l|^ + 2 Iy^)/ 2 (green squares) and l|^ (bine squares) polarizations for H=0.88. l|^ typically probes the magnetic fluctuations 
along c, <M^> whereas (l|^ + 2l|^)/2 is proportional to The background is decreasing with temperatnre as expected 

for a spin incoherent background scattering. Solid lines are guides to the eye. 















Supplementary information 


We here give additional information necessary for the analysis of the polarized neutron data. We 
report polarized neutron measurements on a YBa2Cu306.85 (Tc = 89K) single crystal of mass ~ 9.5 
g. The polarized neutron experiments have been performed on two spectrometers: the triple axis 
spectrometer 4F1 ( Laboratoire Leon Brillouin, Orphee, Saclay) and the multi-detector diffractome¬ 
ter D7 (Institut Lane Langevin, Grenoble). Many details concerning polarized neutron and flipping 
ratio on 4F1 has been already discussed in Refs. in the context of the magnetic intra unit cell 
in high temperature superconducting cuprates. 

For the measurements on 4F1, the incident and hnal neutron wave vector are set to kj = kp = 
2.57A“^ and the energy resolution is ~ 1 meV. On D7, the incident wavelength is taken as Ai = 4.86A 
(= kj = 1.29A“^), which corresponds to Ei = 3.47meV. Since it is a diffractometer, neutrons with 
different hnal energies around Ej are measured. That corresponds to an integration of the magnetic 
huctuations in an energy range between —kpT to Ej, where kp is the Boltzmann constant. 
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FIG. S4: (color online) 4F1 data: Temperature dependencies of normalized magnetic intensity (=Iiuc/Insf where Iiuc is the 
expected magnetic intensity of the lUC order) measured at the wave vector (a) Q = (2,0,0), (b) Q = (0,0,4) as references and 
at (c)-(d) Q = (1,0,L) for the lUC magnetic order study in the P || Q configuration, (c) L = 1, (d) L = 0. The data for Q = 
(1,0,0) (YBa2Cu306.85, full squares) have been calibrated in absolute units using the intensity of the Bragg peak Q = (0,0,4) 
P] and further compared to the YBa 2 Cu 306.6 study (red empty squares) @|. 


Flipping ratio 


The intra-unit-cell (lUC) magnetic order produces a signal at the same position as the nuclear 
Bragg rehections [8|,[^. The magnetic signal is much weaker than the nuclear signal but it acts to hip 
the neutron spin. As polarized neutron measurements enable to distinguish nuclear and magnetic 
scattering, they are a very elegant tool to observe the lUC magnetic order. The scattered intensity 
at a given wave vector Q is then measured in both spin-hip (SF) and non-spin-hip (NSF) channels. 
The hipping ratio, FR = Insf/Isf, dehnes the quality of the polarization. A typical hipping ratio 
is of the order of 40 for 4F1. Furthermore, polarized neutron scattering allows to separate spin 
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orientation along three orthogonal directions. For the conhguration where the polarization P is 
parallel to Q, where the magnetic scattering is exclusively in the SF channel (8|. Consequently, we 
essentially carried out measurements with that polarization on 4F1. 

We show on Fig.EHthe results from the spectrometer 4F1 as a function of the temperature in both 
SF and NSF channels, for the neutron polarization P//Q. To separate the magnetic scattering, it 
is convenient to estimate the inverse of the flipping ratio, which is given by: 


IIFR{T) = Isf/Insf = l/FR^T) + Ruc/Insf (1) 

1/FR° represents the polarization leakage (from the NSF channel into the SF channel), character¬ 
izing the neutron beam polarization quality of the instrument. It is defined at the highest measured 
temperature. 1/FR-1/FR° (=Iiuc/Insf where hue is the expected magnetic intensity for the lUC 
order) represents the normalized magnetic intensity reported in Fig. [SU This quantity is used in 
order to easily compare different studies of the magnetic intensity on different samples. At the hrst 
approximation, 1/FR° is supposed to be temperature independent. However, it appears empirically 
to be slightly temperature dependent [^. This T-dependence is determined by measurements at 
Bragg peaks where the magnetic signal can be ignored such as at large |Q| where any magnetic 
signal becomes vanishingly small. In Fig. [SUa-b, the normalized magnetic intensity obtained at Q 
= (2,0,0) and Q = (0,0,4) is reported where a magnetic signal (if any) is beyond a threshold of 
detection (< 5x10“^). That defines in our experiment, 1/FR‘^(T) as the average of the inverse of 
the flipping ratio at Q = (2,0,0) and Q = (0,0,4) shown in Fig. [SHa-b where no magnetic signal is 
observed, signifying that those two wave-vectors are good reference positions for corrections. 

In order to extract the magnetic intensity, we measured the Bragg reflections Q = (1,0,L) (L either 
0 or 1) as a function of the temperature. As the Bragg intensity appeared on top of a background, 
we also needed to measure the intensity at Q = (0.9,0,L) to get the temperature dependence of the 
background in both SF and NSF channels. Fig. ISSl a shows the temperature dependence of the SF 
and NSF scattering for Q = (1,0,1) and Q = (0.9,0,1). Compared to Bragg reflection, Q = (0.9,0,L) 
is small. That background is then subtracted from the signal at Q = (1,0,L) to keep only the intrinsic 
temperature dependence of the Bragg peak. Fig. [SUb shows the temperature dependence of the SF 
and NSF scattering normalized at the highest measured temperature ~ 260K after the background 
has been removed. The NSF signal has been further divided by the T-dependent 1/FR°(T). One 
can observe an enhancement in the SF channel corresponding to the appearance of the magnetic 
signal below 200K on top of NSF signal normalized by 1/FR°(T). 

Fig. EUc shows normalized magnetic intensity for the Bragg reflections Q = (1,0,1). Starting 
from zero scattering at high temperature, a magnetic signal appears below T^ag ~ 200K ± 20K. 
Our measurements hrst made on the Bragg rehection Q = (1,0,1) are conhrmed by those on Q = 
(1,0,0) (Fig. [SUd) where the data are further compared to a previous study on YBa^CusOfi fiHI. 
Such a value for T^ag agrees with recent accurate determination of T* from resistivity pj. The 
magnetic signal displays a characteristic T dependence that can be htted by (1 — T/T^ag)^^ with 
2/3 = 0.4 ± 0.17 a value similar to what was found in YBa2Cu306.6|ll- 

Once 1/FR— 1/FR° = Iiuc/Insf is determined one can convert hue in absolute units following 
a method described in previous studies[I|, [^. For Q = (0,0,4), the square of the structure factor 
corresponds to ~ 7 barns [l[. All cross-sections are given per formula unit (f.u.). Then, to get the 
Bragg intensity, we have applied the following conversion factor: [lNgF(l,0,L) * 7/Insf(0, 0,4)]. In 
our present case, Insf(1, 0, 0)/Insf(0, 0,4) ~ 4/7, so that, if the amplitude of the magnetic intensity 
for both Q = (1,0,L) is ~ 3 * 10“^ of the nuclear intensity, for Q = (1,0,0) the magnetic intensity 
is of 1.2 mbarn/f.u.. Applying the same procedure for Q = (1,0,1), one gets ~ 3 times weaker 
SF intensity, 0.4 mbarns/f.u.. Compared to previous measurements performed on YBa2Cu306.6 at 
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FIG. S5: (color online) (a) Temperature dependencies of the raw neutron scattering intensities measured at the Bragg peak 
Q = (1,0,1) (full symbols) and at Q = (0.9,0,1) (empty symbols), the background position, in the NSF channel (blue) and in 
the SF channel (black), (b) Temperature dependencies of intrinsic Bragg scattering (background subtracted) at Q = (1,0,1) in 
the SF channel (black) and in the NSF channel (blue) divided by a T-dependent bare flipping ratio FR°(T), so that at high 
temperature there is no magnetic signal. Both pictures shows measurements in the P || Q configuration. 


Q = (1,0,0) 0 and Q = (1,0,1) a. the magnetic intensity in YBa2Cn306.85 is ~ 4 times weaker 
(Fig. [SUd). Following the same hypotheses as before to determine the magnetic moment 0 and 
keeping in mind that the magnetic intensity is proportional to M^, the magnetic moment is here 
estimated to be only ~ O.OS/ib compared to O.lpB for nnder-doped YBCO |l|, 0. 


D7 procedure 


The cnrrent set-np of the cold nentron diffractometer D7 (ILL, Grenoble) has been given in refs. 
nl,[Ii|. It makes nse of a triple-blade monochromator which selects the incoming nentrons. A 


beryllinm filter after the monochromator removes high harmonics. A snpermirror polarizer and the 
Mezei flipper select nentrons with a given spin. D7 is eqnipped with a hxed XYZ polarization mode 
at the sample position. Small gniding magnetic helds are prodnced from a set of electric coils. The 
hnal polarization is then analysed by 66 snpermirror benders over a wide range of scattering angles. 
Two detectors are located after each bender, giving a variable flipping ratio for each detector ranging 
from 15 to 35. 

Both 4F1 and D7 allow the polarization to be arbitrarily oriented. On 4F1, as already described 
in previons papers [3-0, the polarization directions can be adapted to any direction relative to Q. 
Contrary to the 4F1 set-np, the directions of the polarization on D7 are along arbitrary directions 
nnrelated to Q. The polarizations X and Y are in the scattering plane whereas Z is perpendicnlar 
to the plane thns always perpendicnlar to Q. The direction X on D7 is along a direction making 
an angle 24" from ki [121]. Measnrements were performed in both SF and NSF channels, for the 
three X, Y, Z polarizations. We nsed the Lamp software provided by ILL to analyse the data (see 
http://www.ill.en/instrnments-snpport/compnting-for-science/cs-software/all-software/lamp/). We 
applied the standard procednre which consists in a backgronnd snbtraction, a flipping ratio correction 
and a vanadinm normalization to get intensities in absolnte nnits for each channels. 

We carried ont measnrements of the backgronnd coming from an empty sample holder so that 
we conld compare measnrements with and withont the sample. A qnartz rod of the size of onr 
YBa 2 Cn 306.85 sample was pnt in the cryostat in order to get the flipping ratio correction. The 
diffnse scattering from amorphons qnartz is entirely nnclear. Conseqnently, it indicates the ratio 
of the correctly attribnted non-spin-flip scattering cross section to the wrongly attribnted spin-flip 
scattering cross section over all the det^ors (132 detectors). This gives the flipping ratio FR for 
each detector. In the limit of large FR 







11 


jc^r _ ^ 2 ) 

Insf = Insf + ^ (Insf - Isf) (3) 

We apply the flipping ratio correction (as Eqs. |2] and |3]) on our dataset in YBCO using flipping 
ratios determined using the quartz sample. Actually, that correction is insufficient. Indeed, the 
amorphous quartz sample gives broad diffuse scattering and the flipping ratios are generally reliable 
for correcting diffuse scattering. However, systematic errors can occur when measuring different 
scattering, such as an intense Bragg peak which naturally is better collimated and samples a different 
weighted average of the integrated polarization efficiency of the instrument. As a result, we actually 
observed larger flipping ratio on our YBCO single crystal than on the quartz reference. To improve 
the analysis, we made the assumption that no magnetic signal is present at 300K. We then adjusted 
the flipping ratio for each detector so that the scan is featureless at 300K for each Q=(1,0,L) and 
for each polarization. We then applied the flipping ratio correction of Eq. [2] with a corrected set of 
flipping ratios. The intrinsic magnetic signal appears then clearly at lOOK on top of the flat scan 
at 300K (Fig. 2b). Finally, we performed measurements on a vanadium rod, to correct the data for 
detector efficiency and to be able to calculate the cross section in absolute units. Further, following 
the procedure discussed above for 4F1, the conversion in absolute units has been cross-checked using 
the Bragg peak Q=(1,0,0) intensity. 


Neutron cross-sections 


For the case of an ordered state with a magnetic moment pointing along a generic direction M = 
(Ma, Mb, Me), one can estimate the magnetic neutron cross-sections for each neutron polarization Iq, 
with a = (X, Y, Z). For a twinned sample, both directions (1,0) and (0,1) are equivalent and one 

can dehne Mab = y^M^ -|- M^ as the in-plane component and M^ as the out-of-plane component!^. 
The neutron intensity for any arbitrary direction of the polarization a = (X, Y, Z) can be written 
as [Tl[, (^)S = + BGq where BGa is the background. For a XYZ-polarization turned by an 

arbitrary angle (3 between X and Q, can be written as, 

« (1 - gf)M2 -f ~ (4) 

oc + Iz sin^l3 (5) 

oc + Iz cos'll! (6) 

where qi = ^L/|Q| tends to zero in our study for L~ 0. The same formulas can be applied for 
the 4F1 setup but with (3 = 0, the polarization along X for 4F1 being along Q. 

The total magnetic intensity is obtained from the sum of the three cross-section in the SF channel 

_gp _gp _gp 

Ix + ly + Iz whose temperature dependences is reported in Fig.2 of the manuscript for two L. 
One sees that Ix^ -|- — l|^ oc M^^.- Qi —t 0 (as it is the case for L=0.25), Iz oc M^. 
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That enables us to study both components of the magnetic moment M separately as done in the 

_gp 

manuscript. The cross-section in the SF channel essentially corresponds to the out-of-plane 

_gp _gp _gp 

magnetic component, M^, and conversely the quantity + ly ~ Iz directly probes the in-plane 
component Following the dehnition of the magnetic moment, the tilt angle 6, schematically 

drawn in Fig. 3.b, is dehned by the angle of the magnetic moment with respect to the c axis. It is 
related to magnetic components as tan(6*) = 

In the opposite limit of disordered magnetism, one generally considers the cross section for isotropic 
magnetic correlations, i.e. M? = = M^. This corresponds to a model for paramagnetism that is 


extensively applied on D7 [f^ . Imag is written as. 


= 21^^ + 21 f (7) 

The standard D7 full XYZ polarization analysis separates the (para)magnetic intensity indepen¬ 
dently from the background (such as nuclear and spin-incoherent scattering). The only assumption 
is that the background is not polarization dependent. This is readily the case at position away from 
the Bragg peaks. It is worth to note that Eq. [7] determines the total magnetic scattering only for 
paramagnetic scattering. In this limit, it gives Ipara = |g^S(S -|- 1) where g is the Lande factor and 
S is the effective spin. This does not apply for anisotropic magnetic systems. 


Averaging intensities 


In order to get better contrast of the intensities for particular wave-vectors, we averaged intensities 

N 

h from each N detectors using 1=1] h/N. That quantity is actually related to the integrated area 

i=l 

of the peak. The intensity at each detector h {i.e. at a wave vector g*) is a sum of a background 
BGi and a contribution coming from a peak centered at go with a FWHM of Ag and a magnetic 
amplitude Iq. For a Gaussian description of the peak, h is written as 


li = BGi + /oexp(-41n2(gi - qof /Aq'^) 
The averaged intensity T can be written as. 


( 8 ) 


N 


I = J2li/N = BG + Io 


exp{—4:\n2{x — XqY/A q^) 


2=1 




= BG + h 


Ag 


O’ 


qN - qi 


TT 


41n2 


( 9 ) 


(Qn — di) is the width on which we are averaging intensities, BG is the averaged background level 
taken at high temperature near the peak. The quantity plotted in Fig. 3 of the manuscript is T —BG, 
it is proportional to the product of the intensity at the maximum and of the q-width: oc loAq. As 
Ag oc l/^a6; I then typically measures the ratio of the magnetic intenisity to the in-plane correlation 
length: 7 oc lo/^ab- 
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